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Abstract The equilibrium composition of neutron star 
matter is achieved through weak interactions (direct and 
inverse beta decays), which proceed on relatively long 
time scales. If the density of a matter element is per- 
turbed, it will relax to the new chemical equilibrium 
through non-equilibrium reactions, which produce en- 
tropy that is partly released through neutrino emission, 
while a similar fraction heats the matter and is even- 
tually radiated as thermal photons. We examined two 
possible mechanisms causing such density perturbations: 
1) the reduction in centrifugal force caused by spin-down 
(particularly in millisecond pulsars) , leading to rotochem- 
ical heating, and 2) a hypothetical time-variation of the 
gravitational constant, as predicted by some theories of 
gravity and current cosmological models, leading to grav- 
itochemical heating. If only slow weak interactions are al- 
lowed in the neutron star (modified Urea reactions, with 
or without Cooper pairing), rotochemical heating can 
account for the observed ultraviolet emission from the 
closest millisecond pulsar, PSR J0437-4715, which also 
provides a constraint on \dG / dt\ of the same order as the 
best available in the literature. 
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1 Introduction 

Neutron star matter is composed of degenerate fermions 
of various kinds: neutrons (n), protons (p), electrons (e), 
probably muons (/i) and possibly other, more exotic par- 
ticles. (We refer to electrons and muons collectively as 
leptons, I.) Neutrons are stabilized by the presence of 
other, stable fermions that block (through the Pauli ex- 
clusion principle) most of the final states of the beta 
decay reaction n — ► p + I + v. The large chemical poten- 
tials Hi (« Fermi energies) for all particle species i also 
make inverse beta decays, p + I — ► n + v, possible. The 
neutrinos (y) and antineutrinos (&) leave the star with- 
out further interactions, contributing to its cooling (e . g., 
IShaniro fc Teukolskvll983tlYakovlev fc PethicU200l . The 
two reactions mentioned tend to drive the matter into a 
chemical equilibrium state, defined by rj np i = //„ — (x p — 
Hi = 0. 

If a matter element is in some way driven away from 
chemical equilibrium (r) np i ^ 0), free energy is stored, 
which is released by an excess rate of one reaction over 
the other. This energy is partly lost to neutrinos and 
antineutrinos (undetectable at present), and partly used 
to heat the matter. The heat is eventually lost as thermal 
(ultraviolet) photons emitted from the stellar surface. 

The chemical imbalance can be caused by a change 
in the density of the stellar matter. This can in turn 
be prod uced in different ways. The fir st to be consid- 
ered Cbv lFinzilll96.4lFinzi fc WoUlll968ft was stellar pul- 
sation; however, so far no clear evidenc e for this pro- 
cess has been seen. Gravitation al collapse l)Haensell992t 
iGoureoulhon k Haensell Il993|) and mass accretion are 
also possible mechanisms, but in these contexts the non- 
equilibrium heating is probably overwhelmed by the en- 
ergy released through other channels. 

Here we review our work on neutron star heating 
through beta processes in two other contexts, which we 
consider to be the most promising in revealing infor- 
mation about the physics of dense matter and gravi- 
tation. One is rotochemical heating ijReiseneggeil Il995l 
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1997 [lFernandez & R.eiseneggerl200d hereafter FR05;E 



2006) , in which the precisely measurable decrease of tl 



stellar rotation rate, through the related reduction of tl 
centrifugal force, makes the star contract progressivel 
keeping it away from chemical equilibrium. The moi 
speculative gravitochemical heating is based on the h 
pothesis that the gravitational "constant" may in fa< 
vary in time, causin g a similar contra ction or expansic 
of the neutron star ijJofre et al.ll2006l hereafter JRF06 
We refer to our published papers for a detailed discussic 
of ou r methods and formalism (see also lFlores-Tulian fc 
120061) . Here we restrict ourselves to a general, unifie 
discussion of these two processes and their main imp] 
cations. 



2 Time evolution 

The formalism for calculating the evolution of the tem- 
perature and chemical imbalances for the case of roto- 
chemical heating is described in § 2 of FR05. Here, we 
just outline the fundamental equations and the modi- 
fications required in order to treat the gravitochemical 
case as well. The evolution of the internal temperature, 
T, taken to be uniform inside the star, is given by the 
thermal balance equation, 

[L H (Vn P i, T) - L v (rj np i,T) - X 7 (T)] , (1) 



T = 



C{T) 



where C is the total heat capacity of the star, Lh is 
the total power released by the heating mechanism, L v 
the total power emitted as neutrinos, and L 7 the power 
released as thermal photons. Here and in what follows 
(including all figures), all temperatures, chemical imbal- 
ances, stellar radii, and luminosities are "redshifted" to 
the reference frame of a distant observer at rest with 
respect to the star. 

The evolution of the chemical imbalances is given by 

Vnpl = - [AD,l(Vnpe,T) + A M ,l (Vnpe, T)} 

- [BdAVup^, T) + B Mt i(T) npii ,T)] 

— RnplMi2 + C n piG. 



(2) 



The functions A and B quantify the effect of reactions 
towards restoring chemical equilibrium, and thus have 
the same sign of rj np i (FR05). The subscripts D and M 
refer to direct Urea reactions, 



n — > p + I + v, 
p + I — > n + v, 



(3) 



which are possibly forbidden by momentum conserva- 
tion, and modified Urea, 



n + N 
p + l + N 



p + N + I + 77, 
n + N + v. 



(4) 

where an aditional nucleon N m ust participate in or- 
der to conserve momentum (e. g.. lShapiro fc Teukolskvl 




log t [yr] 

Fig. 1 Predicted time-evolution of the surface tempera- 
ture, T s , of a neutron star with rotochemical heating. All 
curves correspond to stars wi th mass M = 1.4 M@ , equa- 
tion of state A18+<5u+UiX* lAkmal et alj lT998l which al- 
lows only modified Urea processes to occur, and magnetic 
dipole spin-down with Bdipoic = 2.5 x 10 11 G. Each curve 
is labeled by the assumed initial rotation period. The up- 
per limits corres pond to observat i onal c onstraints for pul- 
saxs B0950+08 (IZaylin fc Pavlov! l200ll and J0108-1431 
(Mianani et al. 2003:[ as interpreted by IRargaltse v et all 
2004), both of which have magnetic dipole field strengths 
very close to the assumed value. 



Il983t lYakovlev fc Pethickll2004|) . The scalars R np i and 
C np i quantify the departure from equilibrium due to the 
changes in the centrifugal force (oc fifl) and the gravita- 
tional constant (G) , being positive and de pending on the 
stella r model and equation of state (FR05: lReisenegger et alJ 
120061 JRF06). 

Figure^shows the solution of the coupled differential 
equations ^ and [21 for the evolution of a classical pulsar 
with a moderate magnetic field and different assumed 
initial rotation periods under pure rotochemical heating 
{fif} < 0, G = 0). It can be seen that, for very fast 
initial rotation, the pulsar can be kept warm beyond the 
standard cooling time of ~ 10 7 yr, at a level that is close 
to current observational constraints. 

The case of a "millisecond pulsar" (a neutron star 
with fast rotation and a weak magnetic dipole field) is 
illustrated in Figure [3 It first cools down from its high 
birth temperature, while the chemical potential imbal- 
ances r] np i slowly increase due to the decreasing rotation 
rate, until rotochemical heating increases the tempera- 
ture again, and the reactions stop the rise of the chemical 
potential imbalances. 



3 Stationary state 

If the relevant forcing (flfl or G) changes slowly with 
time, the star eventually arrives at a stationary state, 
where the rate at which the equilibrium concentrations 
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log t [yr] log t [yr] 



Fig. 2 (taken from FR05) Evolution of the internal temper- 
ature and chemical imbalances under the rotochemical heat- 
ing effect for a 1.4Mq star calculated with the A18 + 5v + 
UIX* equation of state itAkmal et al. 1998), with initial tem- 
perature T — 10 s K, null initial chemical imbalances, and 
magnetic dipole spin-down with field strength B = 10 8 G 
and initial period Pq — 1 ms. 




log t lyrj 



Fig. 3 (from FR05) Evolution of the internal temperature 
under rotochemical heating for different initial temperatures. 
We set rjnpe = rjnpfj. = 0. The short-dashed line is the quasi- 
equilibrium solution, obtained by solving T = and f)npe = 
"Hnpfi = 0. The stellar model and spin-down parameters are 
the same as in Figure [5] 



are modified by this forcing is the same as that at which 
the reactions drive the system toward the new equilib- 
rium confi guration, with hea ting and cooling balancing 
each other l|Reiseneggerll995fl . The evolution to this state 
for pure rotochemical heating is illustrated in Figure [21 
Figures [3] and 0] show that the state reached is indepen- 
dent on the assumed initial conditions. 

The properties of this stationary state can be ob- 
tained by the simultaneous solution of equations Q and 



Fig. 4 (from FR05) Evolution of the internal temperature 
under rotochemical heating for different initial chemical im- 
balances rj npe = ?7„p M = r) and the same initial temperature 
T = 10 8 K at t = 0. The line styles, the stellar model, and 
the spin-down parameters are the same as in Figure 



P)l with T — rjnpi — 0. The existence of the stationary 
state makes it unnecessary to model the full evolution of 
the temperature and chemical imbalances of the star in 
order to calculate the final temperature, since the sta- 
tionary state is independent of the initial conditions (see 
FR05 for a detailed analysis of the rotochemical heating 
case). For given values of flfl and G, it is thus possible 
to calculate the temperature of an old pulsar that has 
reached the stationary state, without knowing its exact 
age. 

When only modified Urea reactions operate, it is pos- 
sible to solve analytically for the stationary values of the 
photon luminosity and chemical imbalances rj^ l: as 

a function of stellar model and current value of flfl 
and G. The reason for this is that the longer equili- 
bration timescale given by the slower modified Urea re- 
actions yields stationary chemical imbalances satisfying 
Vnpi 3> kT. In this limit, the term Lh — L v in the thermal 
balance equation can be written as KL e r)^ pe + Kl^tj^^, 
where Kl,i are positive constants that depend only on 
stellar mass and equation of state (FR05; JRF06). For 
typical equations of state, the photon luminosity in the 
stationary state is 



/.:' - io ; " 



G/G 



P 3 



3 x 10" 11 



8/7 



erg s 



(5) 



where P ms is the rotation period in milliseconds, and 
P_20 is its time derivative in units of 1CP 20 (dimension- 
less), and the effective surface temperature of the star in 
the stationary state is 



(2 - 3) x 10 5 



P-i 



G/G 



P 3 
ms 



3 x 10 



-li 



yr 



2/7 



K. (6) 
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Finally, the timescale for the system to reach the sta- 
tionary state is 

-6/7 



2 x 10' 



P-5 



G/G 



p3 



3 x 10- 11 yr- 1 



yr. 



(7) 



4 Comparison to observations 

In order to verify this model and constrain the value 
of | G/G |, we need a neutron star that (1) has a mea- 
sured surface temperature (or at least a good enough 
upper limit on the latter), and (2) is confidently known 
to be older than the timescale to reach the stationary 
state. So far, the only object satisfying both conditions 
is the millisecond pulsar closest to the Solar System, 
PSR J0437-4715 (hereafter J0437), whose surface tem- 
perature wa s inferred from an HST- STIS ultraviolet ob- 
servation bv lKarealtsev et aD (120041). Its spin-dow n age, 



5 x 10 9 vr (e.g.. Ivan Straten et al.ll200l[K and 

the cooling age of its white dwarf companio n. t\vd — 
(2.5 - 5.3) x 10 9 yr i|Hansen fc Phinnevlll998h . are much 
longer than the time required to reach the steady state 
for both rotochcmical and gravitochcmical heating, in 
the latter case under the condition that \G/G\ > 10~ 13 
yr" 1 . 

Consequently, we consider stellar models constructed 
from different equations of state, wit h masses satisfying, 
the co nstraint obtained for J0437 bv Ivan Straten et alJ 
l|200ll) . M psr = 1.58 ± 0.18M Q , and calculate the sta- 
tionary temperature for each. Figure[S]compares the pre- 
dictions for the case of pure rotochemical heating with 
the measured spin-down parameters of J0437 (for various 
equations of state and neutron star ma sses) to the tem- 
peratu re inferred from the observation o flKargaltsev et aT] 

In Figure H3 we compare the same observational con- 
straints on the temperature of this pulsar to the theoret- 
ical predictions for pure gravitochemical heating, assum- 
ing \G/G\ = 2 x 10~ 10 yr -1 . As can be seen, this value 
is such that the stationary temperatures of all stellar 
models lie just above the 90 % confidence contour, and 
therefore represents a rather safe and general upper limit. 

When the stellar mass becomes large enough for the 
central pressure to cross the threshold for direct Urea 
reactions, T s drops abruptly, due to the faster relaxation 
towards chemical equilibrium. This occurs in two steps, 
as electron and muon direct Urea processes have differ- 
ent threshold densities (see, e.g., FR05). Conventional 
neutron star cooling models reproduce observed temper- 
atures better when only modified Urea reactions are con- 
sidered (e. g.. lYakovlev fc Pethickll2004UPagell2006j) . Re- 
stricting our sample to the equations of state that allow 
only modified Urea reactions in the mass range consid- 
ered here, namely A18 + 5v, A18 + Sv+ UIX, BPAL21, 
and BPAL31, we obtain a more restrictive upper limit 
on \G\, as shown in Figure[7| yielding \G/G\ < 4 x 10~ 12 
yr" 1 . 




R [km] 



Fig. 5 (taken from FR05) Surface temperature due to roto- 
chemical heating in the stationary state as function of stellar 
radius for diffe rent equations of state, show n as solid lines 
(APR from lAkmal et al.iri99Sl and BPAL from Prakas heFaTI 
1X9811) . for the spin parameters of PSR J0437-4715. Dashed 
lines are 68% and 90% confidence conto urs of the black- 
body fit to the emission from this pulsar llKargaltse v et all 
2004). Bold lines indicate, for each equation of state, the mas s 
range allowed by the constraint of Ivan Straten et alJ (2001), 
Mpsr = 1.58 ± O.18M . BPAL32 and BPAL33 allow direct 
Urea reactions in the observed mass range of PSR J0437. 
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Fig. 6 (from JRF06) Surface temperature due to gravi- 
tochemical heating in the stationary state as function of 
stellar radius for different equations of state. The value of 
\G/G\ = 2 x 10~ 10 yr -1 is chosen so that all stationary tem- 
perature curves lie above the observational constraints. Oth- 
erwise, the meanings of lines and symbols are as in Figure |5] 



5 Discussion and Conclusions 

5.1 Rotochemical heating 

Using the equations of state that allow only modified 
Urea reactions within the allowed mass for PSR J0437- 
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Fig. 7 (from JRF06) Same as Figure [S] but now the value 
of \G/G\ — 4 x 10~ 12 yr -1 is chosen such that only the sta- 
tionary temperature curves with modified Urea reactions are 
above the observational constraints. 



4715, rotochemical heating predicts an effective tempera- 
ture in the narrow range T s>eg = (6.9—7 .9) x 10 4 K, about 
20% l ower than the blackbody fit of Kar galtsev et 

There are three possible reasons why the predic- 
tion does not quite match the observation: 



1 



We are not taking superfluidity into account. This 
would reduce Urea reaction rates, lengthening the 
equilibration timescale and raisi ng the stationary-state 
temperature l)Reiseneggerlll997j) . 
We are neglecting other heating mechanisms (some of 
them directly related to superfluidity), which could 
further raise the temperature at any stage in the ther- 
mal evolution. Nonetheless, in millisecond pulsars, 
all proposed mechanisms appea r to be less impor- 
tant than rotochemical heating <|Schaab et al.lll999t 
IKargaltsev et alJl2004|h 

The thermal spectrum could deviate from a black- 
body, as for the isolated neutron star RX J1856-3754, 
which has a well-determined blackbody X-ray spec- 



to its earlier, faster rotation, which may have built up 
a significant chemical imbalance that is currently being 
decreased by ongoing reactions in its interior (see Fig.^l. 
Depending on its initial rotation period, its surface tem- 
perature may be substantially smaller than both J0437's 
observed temperature and its own current upper limit. 



5.2 Gravitochemical heating 

Table 1 lists some of the many experi ments perfor med 
so far to test the constancy of G (see lUzanl 120031 and 
lwull20oa for recent reviews). The second column con- 
tains the upper limits on its time variation, most usefully 
expressed as \G/G\, and the third is a rough time scale 
over which each experiment is averaging this variation. 
Based on the latter, the experiments can be separated 
into three classes. The first two experiments on the list 
measure the variation of G from the early Universe to 
the present time, and the constraint on the present-day 
value of \G/G\ is based on assuming a time dependence 
G(t) oc t~ a , where t is the time since the Big Bang, and 
a is a constant constrained by these experiments. The 
next four are sensitive to variations over long timescales, 
10 9-10 yr, but without reaching into the very early Uni- 
verse. The last four experiments measure the change of 
G directly over short, "human" timescales of years or 
few decades. Even though results from the first category 
are nominally the most restrictive on a long-term vari- 
ation of G, they depend crucially on the assumed form 
of the variation of G near the Big Bang. Thus, it is still 
useful to consider measurements of the second and third 
categories, which could directly detect variations of G in 
more recent times. 

The new method advocated here, namely gravito- 
chemical heating of neutron stars, falls closest to the sec- 
ond category, as its timescale is much longer than human, 
but does not reach into the early Universe. However, 
it probes somewhat shorter timescales than the other 



WI11U1 llao a vvcii-ucucxmiucu ui&civuuu.y vY-iay opci- mpt.Vinds in this category. In the most general case, when 
trimi that underpredicts the optical flux (Walter fc Mattfe^ ij rca react i ns are allowed to operate 



Il997j) , indicating a more complex spectral shape of its 
thermal emission. 

IKargaltsev et all l|2004[) stress that PSR J0437-4715 
has a higher surface temperature than the upper limit for 
the younger, "classical" pulsar J0108-1431, T. < 8.8 xlO 4 



we obtain 

an upper limit \G/G\ < 2 x 10~ 10 yr -1 . Restricting the 
sample of equations of state to those that allow only 
modified Urea reactions, we obtain a much more restric- 
tive upper limit, \G/G\ < 4 x 10 -12 yr -1 on a time scale 
10 8 yr (the time for the neutron star to reach its quasi- 



K, infe rred from the optical non-detection bv lMTgnani et a l.l stationai T state), competitive with constraints obtained 



(2003) and shown in our Figure [H In the rotochemical 
heating model, these two pulsars are in very different 
regimes: J0437 is in the stationary state in which its 
temperature can be predicted from its spin-down param- 
eters, whereas J0108 has a 680 times smaller spin-down 
power (oc QS7), and will therefore not reach a detectable 
stationary state. Its equilibration timescale, according to 
equation J7J, is 2 x 10 11 yr, longer than the age of the 
Universe and certainly much longer than the spin-down 
age of the pulsar. Thus, its heat content (if any) is due 



from the other methods probing similar timescales. How- 
ever, since the composition of matter above nuclear den- 
sities is uncertain and millisecond pulsars are generally 
expected to be more massive than classical pulsars, we 
cannot rule out the result for the direct Urea regime. 

Further progress in our knowledge of neutron star 
matter will allow this method to become more effective 
at constraining variations in G. The method can also 
be improved with an increased sample of objects with 
measured thermal emission or good upper limits. 
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Table 1 Previous upper bounds on \G/G\. 



Method 


|G/Gj max [1(T 12 yr- 1 ] 


Time scale [yr] 


Reference 


Big Bang Nucleosynthesis 


0.4 


1.4 x 10 1U 


CoDi ct al. (2004) 


Microwave background 


0.7 


1.4 x 10 10 


Nagata et al. (20041 


Globular cluster isochrones 


35 


10 10 


Deel'Innocenti et al. (1996) 


Binary neutron star masses 


2.6 


10 10 


Thorsett (1996) 


Helioseismology 


1.6 


4 x 10 9 


Guenther et al. (1998) 


Paleontology 


20 


4 x 10 9 


Eichendorf fc Reinhardt (1977) 


Lunar laser ranging 


1.3 


24 


Williams et al. (2004) 


Binary pulsar orbits 


9 


8 


KasDi et al. (1994) 


White dwarf oscillations 


250 


25 


Bcnvenuto ct al. (2004) 
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